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ABSTRACT 

High-energy emission (from the X-ray through the 7-ray band) of Flat Spectrum Radio 
Quasars is widely associated with the inverse Compton (IC) scattering of ambient 
photons, produced either by the accretion disk or by the Broad Line Region, by high- 
energy electrons in a relativistic jet. In the modelling of the IC spectrum one usually 
adopts a simple black body approximation for the external radiation field, though 
the real shape is probably more complex. The knowledge of the detailed spectrum 
of the external radiation field would allow to better characterize the soft-medium X- 
ray IC spectrum, which is crucial to address several issues related to the study of 
these sources. Here we present a first step in this direction, calculating the IC spectra 
expected by considering a realistic spectrum for the external radiation energy density 
produced by the BLR, as calculated with the photoionization code CLOUDY. We find 
that, under a wide range of the physical parameters characterizing the BLR clouds, 
the IC spectrum calculated with the black-body approximation reproduces quite well 
the exact spectrum for energies above few keV. In the soft energy band, instead, the IC 
emission calculated using the BLR emission shows a complex shape, with a moderate 
excess with respect to the approximate spectrum, which becomes more important 
for decreasing values of the peak frequency of the photoionizing continuum. We also 
show that the high-energy spectrum shows a marked steepening, due to the energy 
dependence of the scattering cross section, above a characteristic energy of 10-20 GeV, 
quasi independent on the Lorentz factor of the jet. 
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1 INTRODUCTION 

Powerful gamma-ray emission (at energies E > 100 MeV) 
is a distinctive feature of the group of Active Galactic Nu- 
clei collectively known as blazars. As clearly shown by their 
Spectral Energy Distributions (e.g. Fossati et al. 1998), the 
gamma-ray emission belongs to a component different from 
the low-energy synchrotron bump, commonly ascribed to 
relativistic electrons in a jet pointing close to the observer. 
Among the models advanced to account for the high-energy 
component, the most widely discussed (also known as "lep- 
tonic models" to distinguish them from the other models in- 
volving hadronic reactions, e.g. Miicke et al. 2003) assumes 
that the emission is produced through the inverse Comp- 
ton scattering between the synchrotron-emitting electrons 
and either the synchrotron photons (in the Synchrotron Self- 
Compton model, Maraschi et al. 1992) or photons from the 
regions outside the jet (the so-called External Compton mod- 
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els). In turn, external radiation can be dominated by the 
emission of the accretion disk (Dernier & Schlickeiser 1993), 
reaching the jet either directly or through reflection by free 
electrons, or by the emission of the clouds belonging to the 
Broad Line Region (Sikora et al. 1994). The most powerful 
blazars (Flat Spectrum Radio Quasars) display luminous 
broad emission lines and, in some cases, also show the clear 
presence of the "blue bump" , usually associated with the di- 
rect emission from the accretion disk (Sun & Malkan 1989). 
These indications strongly suggest that in FSRQs the high- 
energy component is dominated by the radiation produced 
through the EC process, likely with photons produced by the 
BLR (although the SSC emission can provide an important 
contribution in the medium-soft X-ray band). 

Current models of the EC emission adopts a quite sim- 
plified description of the external radiation field, usually ap- 
proximated with a bumpy, black-body like, spectrum (e.g. 
Inoue & Takahara 1996, Dermer, Sturner & Schlickeiser 
1997, Ghisellini et al. 1998, Bottcher & Bloom 2000). In 
the case of EC models considering BLR radiation, this ap- 



2 Tavecchio & Ghisellini 



proach is justified with the argument that the spectrum of 
the seed photons is dominated by few strong emission lines 
(in particular the Ly a line). However, the actual spectrum 
of the soft photons can be more complex: besides the lines, 
there is the important contribution of the continuum emis- 
sion, which (depending on the parameters characterizing the 
BLR clouds) could imprint important signatures in the EC 
emission. In particular, the shape of the soft radiation field 
would be reflected in the EC spectrum in the medium-soft X- 
ray band, produced by the electrons at the lower end of the 
energy distribution. The determination of the exact spectral 
shape in this band would be crucial for several reasons: 1) 
interesting constraints on the bulk Lorentz factor of the jet 
and the minimum electron energy can be derived from the 
observed spectral hardening below about 1 keV, expected 
to trace the shape of the seed photons (e.g. Tavecchio et al. 
2007; Ghisellini et al. 2007) ; 2) alternatively, if the observed 
hardening below 1 keV in some high-redshift blazars is at- 
tributed to the presence of absorbing material at the quasar 
redshift (e.g. Elvis et al. 1994, Cappi et al. 1997, Reeves et 
al. 1997, Fiore et al. 1998, Fabian et al. 2001a,b, Bassett et 
al. 2004, Yuan et al. 2005, Page et al. 2005), the knowledge 
of the underlying continuum is essential to characterize the 
amount and the physical state of the intervening gas; 3) a 
correct description is also essential in view of the possible 
contribution in this band of the bulk-Comptonized radia- 
tion from cold leptons in the jet (Begelman & Sikora 1987, 
Celotti, Ghisellini & Fabian 2007). 

Apart from these aspects, specifically related to the 
modelling of the soft X-ray emission, the knowledge of the 
diffuse radiation filling the regions surrounding the jet has 
a strong impact in the determination of the opacity of the 
environment to the gamma-ray photons (Liu & Bai 2006, 
Reimer 2007). The assessment of the effects of this "inter- 
nal absorption" is mandatory in view of the possibility to 
use the absorption of gamma-rays to probe the optical-UV 
cosmic background (Chen et al. 2004, Reimer 2007). 

With these motivations, in this work we intend to ex- 
plore some of the implications of the detailed modelling of 
the diffuse radiation field from the BLR for the expected EC 
spectrum of powerful blazars. In particular, we calculate the 
spectrum emitted by the clouds, for a relatively large range 
of the parameters, using the photoionization code CLOUDY 
(Section 2). These spectra are then used to calculate the EC 
spectra (Section 3), which can be compared with the spec- 
tra obtained using the black-body approximation. Finally, 
in Section 4 we discuss the results. 



2 THE EXTERNAL RADIATION FIELD 

2.1 Setting the stage 

The geometry assumed for the modelling is depicted in Fig. 
1. The central accretion flow illuminates the clouds resid- 
ing in the BLR (characterized by the [total: ionized + neu- 
tral] hydrogen density n and the hydrogen column density 
Nh), assumed to be a spherical shell with inner radius R ln 
and thickness AT?. In the calculation we assume that the 
clouds cover a fraction C = fl/^n = 0.1 of the solid angle 
viewed from the central illuminating source. The emission 
from the illuminated face of the clouds (comprising the re- 
flected incident continuum, the emission lines and the diffuse 




Figure 1. Sketch of the geometry assumed in the model (not to 
scale). The (uniform) BLR is assumed to be a spherical shell with 
thickness AR and inner radius i?; n , illuminated by the central 
continuum with luminosity Lyi. See text for details. 



continuum from clouds, see e.g. Korista & Ferland 1998) is 
assumed to isotropically fill the region within i?i n and it is 
calculated with version 05.07 of CLOUDY, described by Fer- 
land et al. (1998)3- For simplicity we discuss only the case of 
solar abundance and we neglect the effect on the clouds of 
the net diffuse emission from the clouds themselves. We as- 
sume that the distance of the jet from the central black hole 
is negligible compared with the radius of the BLR. This as- 
sumption greatly simplifies the geometrical treatment. Fur- 
thermore we do not consider the possible contribution of 
photons originating in the disk, reaching the jet either di- 
rectly or after being re-isotropized by intercloud gas (e.g. 
Celotti et al. 2007). 

We adopt the spectrum of the illuminating continuum 
modelled as a combination of a UV bump with an X-ray 
power- law v usually assumed in this kind of calculations 
(AGN model in CLOUDY, e.g. Korista & Goad 2001 and 
references therein): 



Lm(u) oc v 



exp 



hv 
kTiE 



-Av-\l) 



Above 100 keV we assume a steep {v~ 2 ) power law mimick- 
ing the expected cut-off. This shape reproduces reasonably 
well the average spectra measured from AGNs (Francis et 
al. 1993; Elvis et al. 1994). The bump is assumed to have a 
cut-off in the IR band at 9.1 /im (Korista et al. 1997), since 
most of the IR emission from AGNs is believed to originate 
in warm dust outside the BLR. The normalization constant 
A is fixed to provide a value of a ox — 1.4, as typically ob- 
served in AGNs (e.g. Zamorani et al. 1981). The continuum 
is then fixed by assuming a normalization for the spectrum 
(parametrized using the bolometric luminosity Lm) and the 
temperature of the bump, Tbb- 

Using this model for the continuum we calculate the 
diffuse emission of the BLR for the hydrogen density and 



1 See also http://www.nublado.org/ 



BLR and emission of blazars 3 



the hydrogen column density of the clouds in the ranges 
n = 10 9 - 10 11 cm" 3 and N H = 10 22 - 10 24 cm" 2 , re- 
spectively. These intervals cover the region of the parameter 
space allowed by detailed modelling of the lines observed for 
few selected (radio-quiet) AGNs (e.g. Korista & Goad 2000; 
Kaspi & Netzer 1999). A caveat (which applies also to the 
choice of the spectral shape of the illuminating continuum) 
is that the present knowledge of the geometry and the physi- 
cal state of the gas in the BLR mainly relies on the studies of 
(few) well-observed radio-quiet AGNs (mainly Seyfert galax- 
ies). Indeed it is conceivable that for the sources for which 
we are interested here, namely powerful, radio-loud quasars, 
the conditions could be different. Note moreover that, for 
simplicity, we model the BLR with a unique value of the 
physical parameters, although the recent detailed models of 
emission lines quoted above seem to support the presence of 
a stratified BLR, with the parameters (density, column den- 
sity, covering factor) varying with distance from the central 
source. 

2.2 The spectrum of the BLR 

In Fig. 2 we report some examples of the spectra of the 
diffuse BLR emission [plotted in the vL(y) representation], 
calculated for different values of n and Nh, assuming the 
values Tbb = 10 s K (which fixes the peak of the blue bump 
around ~ 10 15 Hz, dashed lines), Lm — 3 x 10 47 erg/s and 
R in = 10 18 cm. 

The large value of the luminosity is suitable to charac- 
terize the powerful quasars at medium-high redshift. Disk lu- 
minosity L > 10 47 erg/s have been inferred from the optical- 
UV data of sources recently observed with Swift/UVOT (e.g. 
Tavecchio et al. 2007, Sambruna et al. 2007). The radius of 
the BLR has been fixed to a value between those provided 
by the relations found by Kaspi et al. (2007) and Bentz et 
al. (2006). These relations connect the luminosity at a given 
wavelength to the distance of the clouds emitting a given 
line, estimated through the reverberation mapping tech- 
nique. The relations recently derived by Kaspi et al. is based 
on the CIV line and the continuum is evaluated at A = 1350 
A, while the Bentz et al. relation uses the H/3 line and the 
continuum at A = 5100 A. In both cases the relation has 
the same dependence, Rblr oc (ALa) 0,5 . The values of the 
monochromatic luminosities L1350 and L5100 depend on the 
spectral shape of the illuminating continuum. For the shape 
and the luminosity assumed here, the two relations provide 
relatively different values for Rblr,, Rblr — 5 x 10 17 cm 
and Rblr = 3x 10 18 cm for Kaspi et al. and Bentz et al., 
respectively. For the calculation shown in Fig. 2 we assume 
the (logarithmic) average value Ri n = 10 cm. 

It is worth noting that, although calculated for a specific 
set of parameters, the spectra reported in Fig. 2 can be 
easily generalized to different values of Lm and BLR radius. 
Indeed, the normalization of the diffuse spectrum simply 
scales as the luminosity Lm, while the shape (continuum 
and lines) only depends (at the first order) on the ionization 
parameter £ = L a \/nRf n (in terms of £, the spectra in Fig. 
2 correspond to £ = 300, 30 and 3, from bottom to top, 
respectively). Different combination of Lm, i?; n and n can 
thus produce similar spectra. Note that, if Lm and Ri n are 
linked by a relation of the form Ri n oc L^ 5 as discussed 
above, the ionization parameter in different sources would 



only depend on the density of the clouds (assuming similar 
spectra for the photoionizing continuum), since the ratio 
Lm/ii 2 n would be a fix number in different sources. 

As expected, the total luminosity of the diffuse emission 
is dominated by the optical-UV emission, and, in particu- 
lar, by few prominent lines in this band, most notably the 
Ly a lines of H (1216 A) and Hell (303 A). In all cases the 
underlying continuum around the peak (dominated by re- 
combinations and scattering, both Thomson and Rayleigh) 
contributes to few percent of the total luminosity. However, 
the continuum (mainly due to free-free emission) dominates 
below the peak (below ~ 10 15 Hz), making the spectrum 
much broader than the black body approximation. As we 
discuss later, this fact has a direct consequence on the de- 
rived EC spectrum. 

Some trends are clearly visible in Fig. 2. Increasing Nh 
(left to right) has the effect to increase the total emitted 
luminosity and, importantly, to increase the relative con- 
tribution of the optical-UV continuum with respect to the 
emission lines. The increase of the density n (bottom to top), 
instead, has the effect of reducing the ionization parameter £ 
and leads to increase the importance of the recombinations 
(as witnessed by the prominent edges). 

An interesting feature of these spectra is the presence of 
an important X-ray component, originating from the Comp- 
ton reflection on the cloud gas (analogously to the "reflec- 
tion bump" from accretion disks). As expected, the fraction 
of the X-ray continuum (and of the overall continuum) re- 
flected by the BLR increases (from left to right) with the 
column density Nh- Increasing (from bottom to top) the 
density n, instead, has the effect of reducing the reflected X- 
ray luminosity, since the ionization parameter £ (and thus, 
for the same value of the column density, the number of free 
electrons available for scattering) decreases. The decreased 
ionization state of the gas is also responsible for the deepen- 
ing of the depression visible in the soft X-ray band, caused 
by photoelectric absorption. The presence of an important 
diffuse X-ray emission inside the BLR could have important 
consequences for the propagation of 7-rays with E ~ 100 
MeV produced by the jet (Ghisellini et al. 2008, in prepara- 
tion). 

The basic features of the BLR emission do not strongly 
depend on the actual shape of the photoionizing continuum. 
The effect of a different temperature Tbb is illustrated by 
Fig. 3, reporting the SEDs of the BLR reflection calculated 
for the case n = 10 10 cm" 3 , N H = 10 23 cm" 2 and Tbb = 
2.5 x 10 4 K. The overall shape of the SEDs is similar to the 
high-temperature case (Fig. 2, central panel). However, the 
lower temperature has the effect to decrease the relative flux 
of the ionizing photons, leading to decrease the luminosity of 
the emission lines relative to that of the continuum (whose 
peak, moreover, moves to lower frequencies by a factor of 
?s 2). As we show later, the increasing importance of the 
continuum with respect to the lines can significantly affect 
the resulting EC spectrum. 

3 INVERSE COMPTON SPECTRA 

3.1 The external energy density in the jet frame 

To calculate the IC emission produced by relativistic elec- 
trons scattering off the diffuse emission of the BLR we have 
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Figure 2. Examples of the spectral energy distribution of the emission from the Broad Line Region calculated for different values of gas 
density, from n = 10 9 cm _,! (bottom) to n = 10 11 cm~ 3 (top), and column density, from Njj = 10 22 cm~ 2 (left) to Njj = 10 24 cm -2 
(right). The SEDs indicated with dashed lines show the illuminating continuum, for which Tgg = 10 s K. The model assumes a total 
luminosity of the incident continuum of Lni = 3 X 10 47 erg/s and the BLR clouds (with covering factor C = 0.1) arc assumed to reside 
in a spherical shell with thickness AR = 4 X 10 18 and inner radius i?; n = 10 18 cm. 
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T BB =2.5xl0 4 K 
N H = 10 23 cm" 3 
n=10 10 cm" 3 
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Figure 3. SED of the BLR emission (and the corresponding il- 
luminating continuum, dashed line) calculated for a temperature 
of the bump, Tbb = 2.5 X 10 4 K, density n = 10 10 cm -3 , column 
density Nu = 10 23 cm -2 and the other parameters fixed to the 
values used in Fig. 2. Note the more prominent continuum with 
respect to the higher-temperature case shown in Fig. 2. 
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Figure 4. Radiation energy density of the BLR emission (calcu- 



lated for n = 10 cm , column density Njj 



and 



the other parameters fixed to the values used in Fig. 2) trans- 
formed into the jet frame. The bulk Lorentz factor of the jet is 
fixed to r = 15. For comparison we report the corresponding 
energy density in the quasar frame, before the boosting (in this 
case both the frequency and the energy density are measured in 
the quasar rest frame). The dashed line shows the black-body 
spectrum used to calculate the approximate EC emission. 



first to transform the spectrum of the seed photons to the 
jet frame, assumed to move with bulk Lorentz factor F. We 
consider only the photons emitted by the emispheric portion 
of the BLR above the accretion disk, since the disk "occults" 
the other emisphere. 

Because of relativistic aberration, all photons are ob- 
served in the jet frame as coming into a narrow cone of semi- 
aperture ~ 1/r with the axis aligned along the jet speed. 
Within this cone, the radiation field coming from different 
angles are observed to have different frequencies and inten- 
sities, maximum along the axis and minimum at the edge. 
The result is that an external monochromatic radiation field 
with frequency v is observed in the jet frame to have a 
hard spectrum (oc v' 2 , we use the prime to indicate quanti- 
ties evaluated in the jet frame), between v[ — Yv (photons 
coming from the cone edge) and v' 2 — 2Tv (photons at the 
axis of the cone). This can be shown in the following way. 
For an external radiation field with specific intensity 
the corresponding radiation energy density in the jet frame 
is: 

U'(y') = i J l'(u')dQ' = ^ jf " I{v) f^j dn' (2) 

where we use the invariance of the ratio I(u) / v 3 (e.g. Rybicki 
& Lightman 1979) and the integral is extended over the 
angles 0' = arccos [s! in the range — 1/r. Now we change 
the variable of integration from p! to v using: 



Vv' 



r(i -/?/*') 

and we obtain: 
_2tt_ ,. 



Tj3v> 



d/j,' 



dv 
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(3) 



(4) 



where the new limits become v\ — v' /[T(l + /3)] and v-2 = 
v'/T. In the case of a monochromatic external field, I(y) = 
I 5(v — Vo) we obtain: 

2tt/ 



t/V) 



V 2 for Yu < u' < {l + (f)Tv 



(•») 



Vficvl 

One can check that, es expected, U' = J U(v')dv' = 
(7/3) T 2 J U{v)dv ~ T 2 U. 

Using Eq. ((4} we can now calculate the energy density of 
the BLR radiation in the jet frame. An example is reported 
in Fig. 4 for the reference case n = 10 10 cm" 3 , Nh = 10 23 
cm" 2 , Tbb = 10 5 K. The effect of the boosting is clearly 
visible for the Ly a line, which is transformed into a broad 
(extending from v' = TvL Va to u' ~ 2Yvl Vc ,), hard (oc v' 2 ) 
feature. Due to the relativistic broadening most of the lines 
merge in few broad features. In view of the comparison be- 
tween the EC spectrum calculated by using the BLR emis- 
sion and the approximated one, we define the black body 
which best approximate the boosted BLR emission, with 
the peak corresponding to the frequency v' B b = ^-^ v hy a 
and with normalization equal to that of the integral of the 
boosted BLR radiation energy density (dominated by the 
Ly a line). It is shown in Fig. 4 by the dashed line. 



3.2 The EC spectrum 

Using the energy density in the jet frame we can now 
compute the EC spectrum. To this aim we use the radia- 
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Figure 5. EC spectra calculated for the BLR energy densities reported in Fig. 2 {EC/BLR, solid lines) and with the approximating 
black-body (EC/BB, dashed lines). Clearly, the two kind of spectra agree well for frequencies above ~ 10 18 Hz. Below this frequency, 
the EC/BB spectra show a continuum harder than the corresponding EC/BLR spectra. The discrepancy increases for increasing values 
of density and column density of the clouds. 
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tive model fully described in Maraschi & Tavecchio (2003). 
Briefly, we assume that the emission region, moving with 
bulk Lorentz factor V at an angle 9 W = 1/T (therefore the 
relativistic Doppler factoi0 is 8 — T) is isotropically filled 
with relativistic electrons, assumed to have a smoothed bro- 
ken power-law electron energy distribution, extending from 
the energy 7 m i n mc 2 to r y rnax rnc 2 , with slopes n\ and ni be- 
low and above the break located at Lorentz factor 7 P . This 
(purely phenomenological) form has been assumed to repro- 
duce the observed shape of the blazar SEDs, without any 
specific assumption on the acceleration/cooling mechanism 
acting on the particles. 

In Fig. 5 we report the spectra that we obtain for 
all the BLR models shown in Fig. 2 (EC/BLR), using 

7min = l,7p = 100, 7max = 10*, Til = 2, Tl2 = 3.5, T = 15 

(values usually derived from the modelling of FSRQs, see 
e.g. Ghisellini et al. 1998; Maraschi & Tavecchio 2003). The 
electron energy distribution has been normalized arbitrarily. 
For comparison, we also report (dashed line) the spectra cal- 
culated using the black-body approximation defined above 
(EC/BB). This can be simply described by three power laws 
with smooth transitions: 1) an extremely hard power law 
(oc v 2 ) below a "break frequency" Vh t — t^BB^lmin ~ 10 18 
Hz ; 2) a power law with spectral index Qi = (m — l)/2 
between Ub T and v p ~ i^bb^"/ 2 (marking the position of the 
EC peak); 3) a steeper power law with slope a.2 = [n-x — l)/2 
from v p up to the frequency v max . 



3.2.1 The low-energy spectrum 

In all cases the EC/BLR spectra agree quite well with the 
EC/BB ones down to the frequency i-t,r 5 marking the posi- 
tion of the peak of the emission from the electrons at the 
lower end 7 m m. Below this frequency, the slope reflects the 
shape of the spectrum of the seed photons. Thus, the flux 
of the EC/BB decreases as v 2 . Instead, the corresponding 
EC/BLR spectra shows a tail in excess of the EC/BB spec- 
tra, clearly following the excess of soft photons with respect 
to a black-body shown by the BLR continuum (see Fig. 
4). In other words, soft X-ray spectra calculated with the 
EC/BLR tends to be softer than the corresponding EC/BB 
spectra. To better clarify this point we report in Fig. 6 the 
EC/BB (top) and the EC/BLR (bottom) spectra (calcu- 
lated for n = 10 10 cm" 3 and N H = 10 23 cm" 2 ) with the 
contribution of electrons in different "slices" of energy. For 
sufficiently thin slices the EC spectral shape simply reflects 
the incoming spectrum of the seed photons (broadened by 
the boosting into the jet frame). The broad BLR spectrum 
results in a tail (roughly hardening from v° to v 2 ) in excess 
to the narrow peaked black-body shape. 

EC/BLR spectra in Fig. 5 shows clear trends related to 
the features of the corresponding BLR spectra. In particu- 
lar, the luminosity increases for increasing Nh (left to right), 
following the corresponding increase of the luminosity dis- 
played by the BLR emission. Moreover, a better agreement 
between the EC/BLR and EC/BB in the soft X-ray band 
spectra is visible for low Nh- This is clearly related to the 
fact that low column densities imply a minor contribution of 
the continuum to the total BLR luminosity. At large column 
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Figure 6. EC/BB (top panel) and EC/BLR (bottom panel) spec- 
tra (solid lines, for the case n = 10 10 cm" 3 and Nh = 10 23 cm -2 ) 
together with the contribution of electrons in small "slices" of en- 
ergy, 7 =1—2, 2—4, 4-8 and 8—16. For sufficiently small ranges in 
energy these components simply reflect the underlying shape of 
the seed photon distribution. The BLR spectrum, broader than a 
black-body, is the reason of the soft tail displayed by the EC '/BLR 
spectrum. 



densities, instead, the relative importance of the continuum 
grows and therefore the discrepancy between the two kind 
of spectra becomes more and more important. Note that in 
this case even the power law above 1 keV becomes slightly 
softer. 

All the effects described above become more important 
in the case of the low-temperature disk, reported in Fig. 
7. In this C2LS6, cLS already noted, the continuum/line lumi- 
nosity ratio is larger than in the high-temperature case (as 
clearly visible in the spectrum of the boosted BLR radiation 
energy density shown in the inset) and thus the discrepancy 
between the EC/BB and the EC/BLR increases. In partic- 
ular, the EC/BLR continuum appears much broader than 
the corresponding EC/BB spectrum. 



3.2.2 The high-energy spectrum 

At high frequencies both the EC/BLR and the EC/BB spec- 
tra shows a rapid decrease above ~ i o 24-25 Hz. The reason 
for such a cut-off in the spectrum is the decreased rate of the 
IC scatterings in the Klein-Nishina (KN) regime. To derive 
an approximate estimate of the cut-off frequency expected 
from the KN, i/kn, let us assume a monochromatic radi- 
ation field with frequency v . In the rest frame of the jet 
the boosted frequency is v' B ~ 1£v . Electrons with Lorentz 
factors 7 scatter these photons in the Thomson regime as 
long as jhu' B < (3/4:)rn e c 2 , producing photons with fre- 
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Figure 7. EC spectrum calculated for the reference value of the 
BLR parameters and Tbb = 2.5 X 10 4 K. The dashed line reports 
the corresponding EC spectrum calculated with the black-body 
approximation. The insert shows the BLR, radiation energy den- 
sity (in the vU(y) representation) obtained in this case, together 
with the corresponding black body (dashed line). The clear dis- 
crepancy between the EC spectra calculated using the two shapes 
for the radiation energy density is the effect of the minor promi- 
nent "shoulder" displayed by the BLR spectrum below the Ly a 
line. 

quency v' c ~ (3/4)^7 2 . Using this expression for the fre- 
quency of the outcoming photon, the condition that the 
scattering happens in the Thomson limit can be written as 
v' c < 3(m e c 2 /h) 2 lAv' . These photons are observed at the 
frequency vc — 5v' c . Thus, using quantities in the observer 
frame, we can express the condition above as: 



Figure 8. EC/BLR spectra calculated assuming n = 10 10 cm -3 , 
Njj = 10 23 cm -2 , r = S = 15 and, for simplicity, a power- 
law electron energy distribution, extending from 7 m i n = 1 to 
7max = 10 7 and with slope n = 3 (thus producing a spectrum 
with index a = 1) and n = 3.5 (a = 1.25). Solid lines refer to 
the case of the exact Klein-Nishina cross section, while the dashed 
lines report the spectra calculated with the "step" approximation. 
These spectra display a sudden drop of the emission above i^kn 
(dotted vertical line, calculated with u = 2.5 X 10 15 Hz, the 
frequency of the Ly a line) due to the transition into the KN 
regime. In the exact spectra the transition is smoother. 

ellini et al. 1998). In this case a sudden drop of the emis- 
sion, located at the frequency provided by Eq. ©, is visible. 
In the spectra calculated with the exact cross section this 
feature is instead "smoothed" and the spectrum displays a 
continuous decrease. 



35 / m e c 2 \ 2 1 

= 5.6 x lO 24 ^^ Hz = — ^GeV (6) 

Note that typically 5 = T and thus this expression does not 
contain the bulk Lorentz factor of the jet. In this case the 
frequency above which the EC spectrum starts to drop due 
to the KN effects only depends on the frequency of the soft 
photons, v . 

Even if peaked at v , the external photon field is not 
monochromatic, thus, instead of an abrupt break, the EC 
spectrum shows a progressive decrease of the flux. Moreover, 
also the use of the KN cross section contribute to smooth 
the spectrum. An example of the high-energy spectrum (cal- 
culated for the case n = 10 10 cm -3 and Nh = 10 23 cm -2 ) 
is shown in Fig. 8, assuming a simple power-law electron 
energy distribution and two different spectral indices, a — 1 
and a = 1.25 (solid lines). For comparison we also report 
(dashed lines) the same spectra calculated using the "step 
function" approximation for the KN cross section (e.g. Ghis- 



3.3 Application: GB B1428+4217 

As an application of the model discussed above, in Fig. 9 
we show the comparison between the emission calculated 
with the EC model and the XMM-Newton spectrum of GB 
B1428+4217, a powerful blazars located at high redshift 
(z = 4.72). The X-ray spectrum shows a large deficit of 
soft photons with respect to the extrapolation of the power- 
law fitting the data above 1 keV. The lack of soft photons 
in this and other high-z radio-loud quasars has been often 
interpreted as the signature of the presence of large amounts 
of absorbing gas in the quasar environment (Worsley et al. 
2004). However, a direct alternative is that the low-energy 
cut-off is due to the intrinsic hardening of the EC emission 
below f br (e.g. Celotti et al. 2007, Tavecchio et al. 2007). 
Here we tried to reproduce the spectrum with the EC/BLR 
model, trying to find the (probably not unique) set of pa- 
rameters which closely reproduces the data. The correspond- 
ing EC/BB spectrum is also shown for comparison. We fix 
Lji) = 3 x 10 47 erg/s accordingly to that provided by optical- 
UV data, assumed to trace the disk emission. As before, 7?i n 
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Figure 9. XMM-Newton spectrum of GB B1428+4217 (Worsley 
et al. 2004, 2006) (black crosses). The frequency is measured in 
the quasar rest frame, while the observed energy is reported by 
the axis on the top. The red solid line is the EC/BLR model 
reproducing it, using the following parameters: Lm = 10 47 erg/s, 
R BLR = 3xl0 18 cm, n = 10 11 cm" 2 , N H = 10 23 cm" 2 , T = 12.3, 
n\ = 2.2, 7 m i n = 1. The dashed line shows the corresponding 
EC/BB spectrum. 

is fixed close to the average of the values given by Kaspi 
et al. and the Bentz et al. relations. All the other relevant 
parameters have been fixed to "standard" values, n — 10 11 
cm" 2 , N H = 10 23 cm" 2 , F = 12.3, m = 2.2, 7 min = 1. 

As expected, the EC/BLR spectrum shows a tail be- 
low 3 x 10 17 Hz (in the quasars rest frame) softer than v 2 . 
As discussed above, this seems to be an unavoidable and 
robust feature of the EC/BLR spectra, due to the contri- 
bution of the BLR continuum to the EC emission. Instead, 
if the lack of soft photons is the result of absorption, the 
emission should continue to decrease almost exponentially. 
Therefore, observations extending down to very low energy 
and/or sources located at minor distance (so that the spec- 
trum is less redshifted) could be used to to discriminate 
between the alternative interpretation for the deficit of low- 
energy photons. 



4 DISCUSSION 

We have calculated the expected diffuse emission filling the 
BLR region of radio-loud quasars, under quite general as- 
sumptions on the physical parameters of the system. This 
external radiation field is thought to dominate in the jet 
over the synchrotron photons (External Compton emission) , 
hence contributing to the bulk of the inverse Compton emis- 
sion (Sikora et al. 1994). In all cases, few emission lines 
dominate the total luminosity. However, for frequencies be- 
low 10 15 Hz, the emission is dominated by the continuum 
(mainly from free- free emission). We have presented for the 



first time the EC spectra expected from the comptonization 
of realistic BLR spectra, showing that the black body ap- 
proximation (peaked at 1.5 VL Va ) for the external radiation 
field is rather good down to the soft X-ray band. Below 1 
keV, however, the spectra calculated with the BLR spectra 
show an excess, due to the important contribution of the 
BLR continuum. 

The presence of this tail in the EC emission can be 
used to effectively discriminate the origin of the observed 
lack of soft photons observed in several high-redshift radio- 
loud quasars. This evidence is generally attributed to the 
absorption by a huge amount (column densities of the order 
of 10 22 cm" 2 ) of gas in the vicinity of the jet (e.g., Elvis 
et al. 1994, Cappi et al. 1997, Reeves et al. 1997, Fiore et 
al. 1998, Fabian et al. 2001a,b, Bassett et al. 2004, Yuan et 
al. 2005, Page et al. 2005), whose presence is predicted to 
be important especially at high redshifts (e.g. Fabian et al. 
1999). However, the fact that, apparently, only the radio- 
loud quasars display such a deficit of soft photons can sug- 
gest that it can be due to the intrinsic shape of the EC 
emission dominating in this energy band (e.g. Tavecchio et 
al. 2007, Sambruna et al. 2007). Indeed, as we have shown for 
the specific case of GB B1428+4217, this possibility seems 
to be viable even considering a realistic spectrum for the 
BLR emission. However, after the hardening around 1 keV, 
the EC spectrum should display a softer tail extending at 
low energies. Future observations could be effectively reveal 
such a tail. On the contrary, very hard spectra extending 
down into the soft X-ray band can be a strong support to 
the absorption scenario. 

Another interesting point related to the physics of jets 
concerns the value of the minimum Lorentz factor of the 
emitting electrons, 7mm- Having fixed the spectral shape of 
the seed photons for the EC process, the position of the 
break in the spectrum is only related to the product r 2 7^ in 
(see also Maraschi & Tavecchio 2003, Ghisellini et al. 2007). 
Values of 7 m i n larger than 1 would move the position of the 
break at frequencies larger than observed. The break could 
be located at the right position decreasing the value of the 
bulk Lorentz factors. For instance assuming 7 m i n = 2 in the 
case of GB B1428+4217 would require r = 6. These values 
however, seem to be inconsistent with those inferred from 
the modelling of the SED, usually in the range 15-20 (e.g. 
Celotti & Ghisellini 2008). 

The discussion above does not consider the contribution 
to the EC process of the soft radiation coming from other 
regions, as the disk itself (directly or reprocessed) or from 
the putative obscuring torus (Blazejowski et al. 2000). In 
the cases the corresponding EC emissions should extend at 
very low frequencies (close to the optical band), unless the 
value of 7 m in is rather large (e.g. Celotti et al. 2007). There- 
fore in this case the absorption would be the only possible 
explanation for the presence of cut-offs the soft X-ray band. 
Moreover, we also neglect in the discussion the contribution 
of the SSC component, important especially in the soft X- 
ray band. However, from the observational point of view the 
the SSC component seems to be less important in power- 
ful sources (maybe due to an increasing importance of the 
accretion disk and then of the external radiation field: e.g., 
Celotti et al. 2007), leaving an almost "naked" EC compo- 
nent. 

Besides the problems related to the description of the 
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soft EC spectrum, the modelling of the BLR radiation field 
is rather important also in view of the estimate of the opac- 
ity of the quasar environment to 7-rays. Indeed absorption 
of multi GeV photons inside the quasar can represent a seri- 
ous problem for the possible use of blazars as probes of the 
optical-UV cosmic background. Recent works (Liu & Bai 
2006, Reimer 2007) uses approximate spectra of the BLR to 
calculate the expected opacity to gamma-rays. We note here 
that another effect, besides internal absorption, is important 
in shaping the intrinsic high-energy spectrum of powerful 
quasars. Indeed, as we have discussed in Sect. 3.2.2, the EC 
spectrum above few tens of GeV is affected by a strong de- 
pression, due to KN effects. The coupling of these effects 
(KN regime and internal absorption) will probably result in 
a rather soft and dim emission above ft; 10 — 20 GeV, making 
it rather difficult to use blazars as probes of the optical-UV 
extragalactic background, unless some other radiative mech- 
anisms dominates the high energy emission. 
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